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•"*****)  We  have  completed  a  numerical  scheme  for  approximating  the  free 
boundary  of  a  two  dimensional  inviscid  incompressible  jet  impinging 
perpendicularly  upon  a  solid  wall.  A  Gould  plot  shows  good  agreement 
with  an  exact  analytical  method  described  by  Birkhoff  and  Zarantonello 
1957.  Our  method  has  the  advantage  over  the  analytic  procedure  that  it 
allows  for  more  realistic  physical  hypotheses.  An  example  of  a  situation 
to  which  our  method  could  be  extended  is  the  axisymmetric  jet,  for  which 
no  exact  solution  is  know.  It  should  be  emphasized,  however,  that  our 
algorithm  is  an  approximating  scheme;  questions  as  to  the  speed  of  convergence 
and  relative  merits  over  the  finite  elements  scheme  used  by  Sarpkaya  and 
Hiriart  as  reported  by  Gallagher  and  Oden  1975  in  connection  with  the  axi- 
synmetric  jet  remain  to  be  investigated. 

Our  work  should  be  viewed  as  a  first  step  towards  developing  a 
predictive  model  for  the  velocity  and  pressure  distributions  in  the 
impinging  jet  of  a  two-phase  .odium.  The  complexities  of  this  process 
require  alternative  approaches  to  those  followed  in  the  past.  This  search 
for  new  directions  has  been  our  primary  objective.. 
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A  NUMERICAL  SCHEME  FOR  THE  IMPINCilNO  JET 

An  algorithm  is  developed  for  predicting  the  free  boundary 
of  a  2-dimensional  jet  which  impinges  perpendicularly  to  a  solid  wall. 
The  flow  is  assumed  inviscid  and  incompressible. 

The  coordinate  axes  are  chosen  as  in  Fig.  1.  The  ve¬ 


locity  components  are  denoted  u(x,y)  and  v(x,y) .  The  velocity  at  the 

nozzle  exit  is  assumed  of  constant  maanitude  u  .  The  half-width  of 

o' 


Let  ’Mx#y)  be  the  streamf unction .  By  adding  a  constant,  if  necessary, 
we  can  assume  that  i  takes  the  constant  value  zero  on  the  jet  axis 
(which  is  a  streamline).  Vie  then  have: 

v 

( 1 )  Hx,y)  =  /  u(>:,y)Jv 

0 

(2)  |j/(x,y)  -  - v  ( : : , y ) 

[The  second  term  in  the  first  equation  represents  the  flux  across  a 
vertical  line  segment  limited  between  (x,y)  and  (x,0);  we  use  the 
known  fact  that  the  flux  across  <_■  rve  is  the  difference  in  values 


of  the  streamfunction  at  the  endpoints.) 


Algorithm.  It  suffices  to  determine  the  portion  of  the  free  boundary 
that  lies  above  the  jet  axis. 

Choose  a  positive  integer  n.  Set  dx  =  —  ,  and  define  x^  =  dx  , 

xn  =  2dx  ,  etc.  Set  x  =0.  Suppose  the  point  (x  ,y  ) ,  for  some 
t  o  m  m 

m  between  0  and  n,  has  been  found  to  be  on  the  free  boundary,  and 

that  v(x  ,y  )  has  been  determined, 
m  m 


We  first  approximate  ^(x,ym)  by  a  2nd  order  polynomial,  i.c, 


(3) 


^(x,y  )  =  A  (x  -  x  ) “  +  B  (x  -  x  )  +  C 
J  m  m  m  m  m  m 


The  equations  (1)  and  (2)  yield; 


('<  /Y^)  =  u  y 

m  m  go 


^ (1 ,vm)  =  0 

,y  )=  -v(x  ,y J 

ili  Hi  *71  i.i 


C  =  u  y 
m  oJ  o 

=  -71.x  ,y  ) 

‘  =  _  Cm^  5n(n~m)dx 

m  (n-m)2dx2 


Equations  (2)  and  (3)  yield: 


(4) 
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i|»*  (x  .  ,  ,  y  )  =  2  A  dx  +  B  =  -v(x  ,  ,,y  ) 
m+1  Jm  m  m 


m+1 ' 


The  slope  of  the  free  boundary  at  (x  ,y  )  will  be: 

m  J  m 


(5) 


s  = 

m 


v(x  ,y  ) 
m  m 

u  ( x  ,  y  ) 
in  1  ra 


Bernoulli's  equation  on  the  free  boundary  gives  us: 


(O 


u  ( x  ,  y  )  =  /  u  “  -  v  ( x  ,  v  )  1 
m  “m  o  m  *m 


Thus  if  we  know  v(x  ,y  ),  we  can  determine  s  ,  and  approximate 

mm  m 

(7>  Vfl  "  I’m  +  smdx 

[The  above  consists  in  assuming  that  the  free  boundary  is  a  line 

segment:  between  (x  ,y  )  and  {>:  .  ,y  ,  ,  )  .  ] 

m  m  r..  •  -  :n+l 

Vie  make  a  final  approximation,  and  use  (4)  to  get 


(8) 


v  (x 


m+ 1 ' 1 m+ 


)  =  v  ( : 


m+1 


,V  )  =  - (2A  dx 
"  :;i  m 


B  ) 
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and  repeat  the  process.  Note  that  now  wo  have  a  now  point  (x  ,  .  ,  v 

1  L  m+L  'm+1 

on  the  free  boundary,  and  the  vaLue  v(x  . ,v  .). 

m+1  *  m  «■  1 

Since  the  initial  values  x  ,  y  ,  v(x  ,v  )  =  0  are  known,  we  can 

o  o  o-o 

then  obtain  a  sequence  of  points  on  the  free  boundary:  (x^,y^), 
(X2#Y2)»  etc.  The  accuracy  of  our  approximations. 


as  can  be  noted, 


improves  when  dx  is  taken  smaller  (or  equivalently,  when  n  in  large). 
The  process  will  stop  when  v  is  larger  than  uq,  for  then  (6)  is 
undefined. 

Comments .  The  above  should  be  viewed  as  a  first  step  toward?  dcvelopin 
a  predictive  algorithm  for  the  velocity  and  pressure  distributions  in 
the  jet  flow.  A  requirement  we  have  imposed  on  our  model  is  that 
it  be  sufficiently  simple  and  flexible  to  allow  additional  physical 
considerations  (wall  friction,  etc.). 

A  Gould  plot  shows  good  correlation  with  the  hodographic  method. 
Relative  merits  of  our  scheme  over  Finite  Elements,  or  other  approx¬ 
imating  schemes,  remains  to  be  inve stigated . 
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